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^ ABSTRACT 

\^ \ Recent observations suggest the existence of a white dwarf population in the 

H \ Galactic halo, while others suggest that deuterium has been astrated in systems 

at high redshift and low metallicity. We propose that these observations could 
be signatures of an early population of intermediate-mass stars. Such a popu- 
lation requires a Population III initial mass function different from that of the 
solar neighborhood, as perhaps also suggested by the observed cosmic infrared 
background. Also, to avoid overproduction of C and N, it is required that the 
Z = yields of these stars have low (~ 10~^ solar) abundances as suggested by 
some recent calculations. Under these assumptions, we present a model which 
reproduces the observed D vs Si trend, and predicts a high cosmic Type la su- 
pernova rate, while producing a white dwarf population that accounts for only 
~ 1.5% of the dark halo. This scenario can be tested by observations of the 
cosmic supernova rate, and by confirmation and further studies of the putative 
white dwarf halo population. 
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Introduction 



One of the key parameters in the standard big bang model is the density of ordinary 
baryonic matter, expressed as a fraction of the critical density, Qb (often quoted in combina- 
tion with the Hubble parameter as Vtsh'^)- Big bang nucleosynthesis (BBN) has traditionally 
been the primary theoretical tool for determining Qsh"^, and has recently been complemented 
by detailed measurements of the cosmic microwave background anisotropy. These results are 
to be compared to direct surveys of cosmic baryonic matter. While baryons are readily ob- 
servable in stars or stellar remnants, galaxies. X-ray-emitting hot gas, as well as intergalactic 
gas, the low-redshift baryonic budget (Fukugita, Hogan, & Peebles (1998)) is uncertain, and 
leaves some room for dark baryonic objects. Conversely, at higher redshifts, {z = 2 to 3) 
a new estimate of fi^/i^ based on the measurement of Lyman alpha forest intergalactic H I 
absorption lines (and also He II lines) in the spectra of remote quasars has been done and is 
in good agreement with Jl^/i^ derived from BBN (in the range fi^/i^ = 0.01 to 0.03, Riedeger 
et al. (1998), Wadsley et al. (2000)). This is between 3 and 20 times the observed baryon 
density in the luminous parts of galaxies and their presumed precursors, the damped Lyman 
alpha intergalactic neutral hydrogen clouds, amounting to Qsh'^ = 0.0015 to 0.003 (Peroux 
et al. 2001). 

At the local scale, MACHOS, with masses of about 0.5 Mq and detected by gravitational 
microlensing towards the LMC, are possible candidates for some of the dark baryons, and 
amount to ^ 20% of the dark halo (Alcock et al. 2000). Recently, old halo white dwarfs 
may have been detected in sufficient numbers to amount to ^ 2% of the dark halo. These 
white dwarfs, if halo members, represent a large increment, by a factor ~ 10, relative to disk 
white dwarfs, and hence would have a profound impact on Galactic chemical evolution. In 
this paper we discuss the implications of the discovery of dark baryons in the form of old 
halo white dwarfs for BBN, and in particular for the deuterium abundance. 

The standard model of BBN has enjoyed great success in its ability to predict the primor- 
dial abundances of the light element isotopes, D, ^He, "^He and ^Li. Indeed, in the standard 
model, there remains only one remaining unknown parameter, the baryon-to-photon ratio, 
V = UB/n^, which can be expressed as rjiQ = lO^^rj = 274^2^/^^. Roughly, for a value of r]io 
between 1.7 and 7.4, (corresponding to ^2^/^^ between 0.006 and 0.027), we obtain concor- 
dance between the BBN predictions and the observational determination of the abundances 
(Nollett & Buries 2000; Vangioni-Flam, Coc, & Casse 2001; Cyburt, Fields, & Olive 2001). 
In detail, however, there are small but perhaps not insignificant differences in the determi- 
nation of Qsh'^ from the different isotopes. ^He and ^Li concordance is typically achieved 
for lower rjio whereas D/H is best ffi at higher values. To be sure, not every isotope serves 
as an equally good baryo meter. In fact, because of the weak logarithmic dependence on ?7io. 
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■^He is the worst element to use to fix the baryon density. While the ^Li abundance strongly 
depends on rjio, the theoretical uncertainties which stem from uncertainties in the input nu- 
clear cross-section make it difficult to derive precise values of the baryon density. In addition, 
because of the functional dependence of ^Li on rjio, a given observational abundance of ^Li 
gives us two predictions for tjiq. D/H, on the other hand, is both strongly dependent on rjio 
and is relatively free from theoretical uncertainties. One should bear in mind however, that 
there are far more observational determinations of ^He and ^Li at the present time than there 
are of D/H from quasar absorption systems. Nevertheless it is the range in D/H, reported to 
span up to a factor of 20, that has been the center of much controversy over the past decade 
because of the strong sensitivity to the baryon density. The reported D/H abundance has 
now converged towards the lower end of this range, but recent measurements reveal a spread 
by up to a factor of 3 that may reflect intrinsic variations in the D/H abundance. These vari- 
ations would most plausibly arise from stellar processing of primordial gas, or possibly from 
non thermal processes such as photodisintegration by gamma rays, or spallation reactions 
induced by energetic particles. 

The recent results BOOMERanG yield Qsh'^ = 0.02llg:[j[|^ (Netterfield et al. 2001), in 
very good agreement with the estimate from D/H. This is also in agreement with the recent 
results of DASI (Pryke et al. 2001), who find VLbH'^ = 0.022+°;°°^ 

There have also been several new observations of D/H from quasar absorption systems. 
In addition to the earlier determinations by Buries and Tytler (1998a,b) in two Lyman limit 
systems of 4.0 ± .65 and (3.25 ± 0.3) x 10~^ and the upper limit in a third Lyman limit 
system by Kirkman et al. (1999) of 6.8 xlO^^, there have been three new measurements of 
D/H in Damped Lyman a systems by O'Meara et al. (2001): 2.5 ±0.25, by D'Odorico et al. 
(2001): 2.25 ± 0.65, and by Pettini and Bowen (2001): 1.65 ± 0.35 all times 10"^ This last 
value is uncertain due to the poorly defined level of the continuum and the possible blending 
of lines (P. Petitjean, private communication). We note that there is also a measurement 
of D/H in a relatively low redshift system with significantly higher D/H (Webb et al. 1997; 
Tytler et al. 1999). This system has recently been criticized (Kirkman et al. 2001) and we 
will not consider it further here. 

Is there a real dispersion in D/H in these high redshift systems? We note two possible 
trends that go in a similar direction. The data may show an inverse correlation of D/H 
abundance with Si. This may be an artifact of poorly determined Si abundances, or (as 
yet unknown) systematics affecting the D/H determination in high-column density (damped 
Lyman-a, hereafter DLA) or low-column density (Lyman limit systems) absorbers. On the 
other hand, if the dispersion in D/H is real, it has profound consequences, as it must indicate 
that some processing of D/H must have occurred even at high redshift. The second trend is 
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that the data may show an inverse correlation of D/H abundance with HI column density. 
If real, this would suggest that in the high column density DLA systems, which are most 
likely to have undergone some star formation, some processing of D/H must similarly have 
occurred, likewise at high redshift. 

It is interesting to note that local D values measured in the interstellar medium also 
appear to show a considerable dispersion. Linsky (1998) summarized the observations then 
available and argued for a constant value in the local ISM of (1.5 ± 0.1) x 10~^. Recent 
data (see Vidal-Madjar 2001, for a review) are consistent with this mean value, but suggest 
significant variations of the D/H ratio spanning the range 5 x 10~^ to 4 x 10~^. Thus, if D 
shows such variability, perhaps correlated to some unexpected physical processes leading to 
its destruction, one can, in turn, question the abihty of cosmological clouds to represent the 
genuine primordial abundance. 

Thus, a simple average of D/H abundance determinations does not make sense, at least 
without a proper enlargement of the error in the mean due to the poor x^ that such a 
mean would produce. Moreover, if deuterium destruction has occurred, we must question 
the extent to which any of these systems determine the value of VLb- 

There is one obvious objection to the hypothesis that nucleosynthetic processing might 
have occurred in the high-redshift systems. The absorbing systems have low metallicity. We 
now point to seemingly unrelated observations that indicate that primordial processing may 
have occurred, and with substantially less enrichment than is associated with standard yields 
for a solar neighborhood initial mass function (IMF). 

A substantial population of low luminosity white dwarfs has recently been discovered 
(Oppenheimer et al. 2001; Ibata et al. 2000). These objects were identified by their large 
proper motion, indicating halo kinematics. Oppenheimer et al. find that if the objects are a 
halo population of white dwarfs, then they represent a fraction > 2.5% of the dark halo mass. 
There has been some question about whether these objects do reside in the halo (Hansen 
2001; Reid, Sahu and Hawley 2001), but a recent re-analysis of the data by Koopmans & 
Blandford (2001) also concludes that these are halo objects, and finds a similar mass budget 
to that of Oppenheimer et al. We now address the potential implications of any observed 
dispersion in D/H in high redshift quasar absorption systems for a relationship between the 
possibility of D/H destruction in these systems by intermediate mass stars and the ensuing 
production of a population of halo white dwarfs. 

We summarize the importance of D as a baryometer, and its relation to the other light 
element abundances, in §2. In §3 we argue that the observed dispersion in high-redshift D/H, 
if real, cannot be accounted for by standard chemical evolution models. Thus, an explanation 
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of such a dispersion will necessarily involve nonstandard scenarios, involving, e.g., variations 
in the IMF, or outflows of material. We search for models which can explain this dispersion 
as well as the white dwarf observations. We find that such models can be constructed, 
assuming an initial burst of intermediate-mass star formation. It is also necessary that the 
C and N yields of zero metallicity stars are small, as suggested by some recent models. The 
prospects for testing this model and the implications for BBN are summarized in §4. 



2. Determinations of VIb 

As discussed above, each of the light elements can be individually used to determine 
fiBh^ with varying accuracy. For example, a helium mass fraction of Y = 0.238 corresponds 
to rjio = 2.4 or Qsh'^ = 0.009. However, the systematic uncertainties in the ^He abundance 
make it difficult to exclude values of Y as high as 0.25 (see for example, Olive & Skillman 
2001) which because of the weak dependence on r] corresponds to tjiq = 7.4 or Qsh"^ = 0.027. 
With ^Li one can do much better. A plateau value of ^Li/H = 1.2 xlO^^*^ corresponds to 
rjiQ = 2.4 and 3.1. However, the slightly higher value of ^Li/H = 2 xlO~^° (consistent within 
systematic uncertainties) yields riio = 1.7 and 4.5. While there are occasional attempts to 
further increase the '''Li abundance on the basis that some '''Li has been depleted (Vauclair & 
Charbonnel 1998; Pinsonneault et al. 1998,2001; Theado & Vauclair 2001), it is very hard to 
do so, due to recent ^Li observations which leave little room, if any, for ^Li depletion (Fields 
& Olive 1999; Vangioni-Flam et al. 1999) and the lack of dispersion in the Li data (Bonifacio 
& Molaro 1997; Ryan, Norris, & Beers 1999; Ryan et al. 2000). 

In contrast, the strong dependence of D/H on t] makes deuterium a very good baryome- 
ter, provided we have an accurate determination of D/H. The values of D/H measured in 
the high redshift systems however cover the range ?7io = 4.8 for D/H = 4 xlO~^ to r/io = 8.5 
for D/H = 1.65 xlO"^. Thus if the dispersion in the D/H is real (beyond observational 
error as the quoted error bars imply), then one can question the extent that any of these 
determinations truly represent the primordial D/H abundance. 

In the event that the dispersion is real, some amount of D/H must have been destroyed, 
presumably by stellar processing. Since D/H(t) is monotonically decreasing in time, the 
primordial value should be greater than or equal to the maximum of the observed values. 

Stellar processing must involve short-lived stars of intermediate or even larger mass. 
The intermediate mass stars are the proposed white dwarf precursors. However, even an 
IMF peaked near 2-8 Mq could include some massive stars above 10 Mq that would explode 
as SNII and produce the heavy element enriched observed in the absorption systems. An 
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inverse correlation of D with Si abundance is a natural consequence. Since these enriched 
systems are preferentially the more massive and advanced with respect to disk formation, 
as we discuss below, they might also have higher gas column densities and account for the 
additional inverse correlation of D with HI column. 



3. Destruction of D/H by galactic chemical evolution 

The high column density quasar absorption systems are generally considered to be pro- 
togalactic disks. In particular DLAs have the kinematics and ionization characteristic of 
the interstellar medium of massive protodisks. Lyman limit systems have a greater extent 
and display higher ionization than the DLAs. Star formation with a normal IMF could not 
possibly have occurred in these systems to any substantial extent. Typical metallicities are 
0.1 to 0.01 solar, and this is consistent with these systems being predominantly gaseous. 
However in order for there to be enough halo white dwarfs to account for the recent ob- 
servations, the primordial IMF must have been very different from the present day IMF, 
and been dominated by intermediate mass stars. These stars would have formed during the 
protogalactic phase, and could have destroyed significant amounts of primordial deuterium. 

Typical stellar D/H depletion mechanisms usually predict the production of heavy ele- 
ments in closed box models. Indeed, by the time the metallicity reaches 0.1 solar abundance 
in a standard galactic evolution model, only a few percent of the initial D/H ratio has been 
affected. This is the main reason one is led to associate the observed D/H ratio to that of 
the primordial deuterium abundance. It is conceivable, however, that an early generation of 
stars was formed with masses only in the intermediate mass range of 2-8 Mq. These stars 
would not be significant producers of heavy elements (oxygen and above) usually ascribed 
to Type II supernovae. Motivation for this suggestion comes from the discovery by Oppen- 
heimer et al. (2001) of a population of old white dwarfs in the halo of the Galaxy (see also 
Ibata et al. etal 2000; Chabrier 1999; Gates & Gyuk 2001). 

Previous studies of white dwarf halos have examined the consequences when the entire 
halo consists of white dwarfs. There are many objections to this extreme hypothesis, ranging 
from overproduction of He, excessive halo Type la supernova rates, excessive light element 
production in the early (proto-) Galaxy, and excessive C and N production. Constraints 
on the high-redshift population, and possible observable signatures have been studied by 
several groups (Ryu, Olive, & Silk 1990; Fields, Freese & Graff 2000; and references therein). 
The principal refinement we make to these ideas is to note that the old white dwarf halo 
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population is sub-dominant, possibly amounting to as little as 2 percent of the dark lialo.^ 
Nevertheless, the effects on galactic chemical evolution are profound. 

To see this, let consider the mass budget for deuterium destruction. In order to achieve 
noticeable D depletion, one must process a significant amount of material though an early 
generation of stars. A typical galaxy may have Mb = I^^^Mq of visible baryons, mostly in 
the form of stars, and about 10^'^Mq total, the balance being dark matter. We point out 
that a white dwarf halo mass fraction of even 2% can have significant consequences for D 
destruction in protodisk and disk gas. This yields about 2 x lO^^M© in white dwarfs, and 
hence ~ lO^^M©, comparable to the present disk mass, in precursor stars. All 8 x lO^^M© 
of the ejecta is totally astrated. Hence present-day astration factors of 2 or more are easily 
achievable with a white dwarf halo mass fraction of 1% or less, in an interstellar medium 
containing some 6 x IO^Mq in atomic and molecular gas. 

The requirements of D astration and white dwarf production, all with little heavy el- 
ement production, lead to very strong constraints on the required IMF. To avoid heavy 
element production through Type II supernovae, we must demand that the IMF is restricted 
to the 2-8 Mq range (Adams & Laughlin (1996)). At the m < 2Mq end, the constraint is 
that too many stars would be luminous today or would have been seen as a huge population 
of red giants at z ~ 3, as t ^ 2 Gyr at this redshift in an i7m = 0.3 and Vt\ = 0.7 universe. On 
the m > 8Mq end, the problem is that we do not want to make a lot of metals (oxygen and 
up) from supernovae. This constraint is not as strong, and we need to add a few supernovae 
to get some metal production for the absorption systems. 

A new result which independently demands an altered IMF at high redshift comes from 
the far IR background. Models of the FIRB (by ULIGs at high z) overproduce K light today 
(ie too many low mass stars) if a standard IMF is adopted (Cole et al. 2001). An IMF 
peaked at intermediate masses offers an attractive solution: it avoids the low mass stars 
and overproduction of heavy elements but gives the light. The FIRB requires vi^ to be 
about equal to what is produced with a standard IMF for the optical/near-IR background, 
i.e., about 10 — 15 nW/m^ (and to be produced in dust-shrouded ULIGS; see, e.g., Hauser 
and Dwek 2001 for a recent review). The white dwarfs that the ULIGS (which are known 
to be merger- induced bursts of star formation) produce would then (possibly after another 
merger or if they formed in outfiow-boosted star formation) be the halo white dwarfs. The 
kinematics are a weak point but the moral is striking: a peaked IMF may be needed for 



^It has been argued (Reid et al. 2001)) that the detected white dwarfs might be a thick disk population. 
If this is the case, the total mass of white dwarfs is difficult to determine precisely but is considerably smaller 
and would thus not have important consequences for deuterium. 



completely different reasons. 

The obvious objection is that the hypothesised precursor stars will overproduce C and 
possibly N that would be observable in Galactic or extragalactic material (e.g., Fields, Freese, 
& Graff 2000). This is a serious issue due to the nature of the observed C and N abundances 
in primitive systems. For our purposes, the most constraining C and N values are those in 
the same QSO absorbers for which D has been detected. Namely, C has been detected in 
the best-measured Lyman limit systems (Buries & Tytler (1998a; 1998b); O'Meara et al. 

(2001)), and in all three cases [C/H] 2, and [C/0] < or [C/Si] < 0. Results for N are 

similar. Thus, the observations apparently demand that no significant C or N enrichment 
(over that of a-elements) takes place. 

These observations place key constraints on a putative intermediate-mass Population 
111. If the yields of zero-metallicity intermediate mass stars are comparable to those of their 
low-metallicity counterparts (e.g., van den Hoek & Groenewegen (1997)), then indeed, the 
resulting enrichment is extremely large if most of the baryons are processed through white 
dwarfs. Although the enrichment is reduced when only a fraction of all baryons are processed, 
even here, to make a white dwarf halo at 10% of the total leads to C and N abundances 
at about [C, N/H] = —1. As this violates the observational constraints, we conclude that a 
halo white dwarf population cannot be created at high redshift ii Z = intermediate-mass 
stars create significant C and N. It then also follows that the dispersion in D due to chemical 
evolution is small; a similar conclusion was reached by Prantzos & Ishimaru (2001). 

It is far from clear, however, that zero-metallicity intermediate-mass stars behave in this 
way. In fact, the nucleosynthetic yields for these stars are quite uncertain. In this mass range, 
at all metallicities, the yields depend sensitively on the details of several "dredge-up" events 
in which convective mixing brings nuclear burning products up to the envelope of the star. 
These have been studied in detail for solar metallicity stars and for moderately metal-poor 
stars, and the main results of different groups are in reasonable qualitative and quantitative 
agreement. However, there remain basic qualitative uncertainties in the post-main sequence 
evolution of zero-metallicity intermediate-mass stars. 

Specifically, it is as yet unclear — i.e., model-dependent — whether the last and most im- 
portant of these mixing events (the "third dredge-up" ) occurs. There has been considerable 
effort recently to answer this question, but the studies have given conflicting results. Chieffi, 
Dominguez, Limongi, & Straniero (2001) do flnd such a phase, and conclude that zero- 
metallicity intermediate-mass stars should be important sources of C and N. The work of 
Marigo, Girardi, Chiosi, & Wood (2001) does not flnd a third dredge-up, but the authors 
are uncertain whether such a phase would have occurred had they run their models longer; 
what is clear is that the C and N envelope abundances prior to this phase are small (< 10~^ 
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of solar) and so would not create a significantly signature. Finally, Fujimoto, Ikeda, & Iben 
(2000) argue that for stars with M > 4M0 and [Fe/H] < —5, there is no third-dredge up at 
all. 

It thus appears that the question of C and N abundances in Z = stars is both complex 
and unsettled. To simultaneously explain a halo white dwarf population and D depletion at 
high redshift requires that these stars do not eject significant C and N. For the moment, this 
amounts to an assumption, but one with some support from ongoing studies. 

To make the proposed scenario concrete and quantitative, we have constructed chemical 
evolution models which lead to significant D astration at low metallicity and also have 
significant white dwarf production. We will present one such model; of course, given the 
freedom to adjust the IMF (and star formation history), a large class of models are allowed. 
However even with this freedom, the acceptable models share common traits and make 
similar predictions. We have explored the case of a closed-box model. The stellar creation 
or birth function C{m, t) encodes the rate and mass spectrum of new star birth via 

dN^ = C{m,t) dmdt (1) 

We adopt the "bimodal" form C{m,t) = ■?/'i(t)0i(m) + '?/'2(^)02(^^)- We choose the star for- 
mation rate tpi = (Mi/ri)e^*/^^ to be a short-lived burst of star formation with ri = 100 Myr 
and Ml = 0.75M^, with an associated IMF which is sharply peaked at intermediate mass 
stars; this is the population which destroys deuterium and creates white dwarfs. Following 
Adams & Laughlin (1996), we take the IMF to be lognormal 

0i(m) ocm~^exp[— In (m/mc)/2a'^] (2) 

The model we present here has a characteristic mass scale rric = 5Mq and a dimensionless 
width 0" = 0.07 (this is very similar to a gaussian with mean n = 5Mq and width ctg = crmc = 
O.SSMq). Other parameter choices are possible; for higher rric, the needed a is narrower. 
The other term in the creation function describes standard, quiescent star formation. We 
have ^2 = A(7(t)Mgas, where A = 0.3 Gyr~^, and the factor g{t) = 1 — e"*/^-^'^^'' allows for a 
smooth transition from the bursting to the quiescent star formation rate. The IMF 02 is a 
power law with slope 1.35 (Salpeter). 

Results from the model appear in Figure 1. In panel (a) we see the evolution of D/H 
versus metallicity [0/H] (as O is an a-element which is co-produced with Si, we take [0/H] = 
[Si/H]). Three models are plotted. The solid curve corresponds to a bimodal model in which 
we have taken a primordial value D/H = 5.0 x 10"^. The two phases of star formation are 
clear, as each leads to a characteristic drop in D, the burst phase acting at low metallicity and 
the quiescent phase acting at high metallicity. The points are the observed data discussed 
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Fig. 1. — Results from the bimodal star formation model, (a) D evolution as a function 
of metallicity. Data and model curves are as noted in the text; the protosolar value is that 
of Geiss & Gloeckler (1998); the ISM value and errorbar of Linsky (1998) appears as the 
solid line, while the range reported by Vidal-Madjar appears as the dotted errorbar. (b) D 
evolution as a function of remnant production as measured by the total Galactic mass M^^m 
in remnants, almost all of which are white dwarfs, (c) The predicted cosmic supernova rates, 
in arbitrary but comoving units, as a function of reshift. 



earlier. The values from quasar absorption line systems have an arrow indicating the Si 
abundances are lower limits.^ The dashed curve in Figure 1 corresponds to a bimodal model 
with D/H = 4.0 X 10^^. In this case, there is less astration needed at low metallicity, so the 
amount of gas processed in the Pop III burst is lower, i.e.. Mi = O.SM^ (as is the resulting 
mass of halo white dwarfs). Note that the two model curves join once the D has been 
destroyed in the Pop III phase, and follow the same quiescent evolution afterwards. We see 
that in both cases, the curves provide a reasonable fit to the data. Therefore, it is difficult 
to use the existing D/H data and chemical evolution to accurately predict the primordial 



^This arises because the Si abundances are measured in the gas phase, and do not include any Si which 
has been depleted onto dust grains. The presence of grain depletion complicates attempts to interpret 
measurements and trends in high- redshift abundances, and to make matters worse, the depletion levels 
found in the Milky Way (more than an order of magnitude) are apparently not always appropriate for high- 
rcdshift systems. As shown for the case of DLAs (e.g., Pettini, Ellison, Stidel, Shapley, & Bowcn 2000; 
Molaro et al. 2000; Prochaska & Wolfe 1999; Lauroesch, Truran, Welty, & York) Si/Zn can be much higher 
than in the ISM. Si/Zn is typically about 1/3 (i.e., 0.5 dex), and varies from system to system, sometimes 
exceeding the solar ratio. 
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D/H value and hence Qsh'^- Finally, the dot-dashed curve corresponds to a model in which 
there is no Pop III and a primordial value D/H = 3.0 x 10~^ as in O'Meara et al. (2001). In 
this "standard" case there is essentially no D depletion at low metallicity and, as we have 
noted, the dispersion in D is unexplained. 

The deaths of the Pop III intermediate-mass stars not only release D-free material, but 
also produce white dwarfs. This correlation is illustrated in Figure 1(b); we plot the D 
abundance versus the total mass in compact remnants, essentially all of which is in white 
dwarfs. The model used here is that shown in the solid curve of Figure 1(a). We see that the 
burst phase, during which D/H drops to about 2.1 x 10~^, puts about 15% of the baryons 
(about 1.5 X IO^^Mq) into white dwarfs. The quiescent phase adds another IO^^Mq. Note 
well the distinction between this work and previous work. Previous work was centered on the 
possibility of a white dwarf halo. Here we are interested in D/H and producing a few white 
dwarfs. A halo population of a few percent would dramatically ease the previous constraints. 
Indeed, Figure 1(b) shows that the 2.5% observation of Oppenheimer is more than sufficient. 

Note that by showing the extragalactic D vs. Si data on the same plot as the Galactic 
points, we implicitly assume that the Galaxy is a system representative of the quasar ab- 
sorbers. It then would also follow that the white dwarf and carbon production also would 
apply to these systems. However, one need not take the model so literally as to demand 
that all baryonic matter-or even all the baryons in protogalaxies-foUow this evolution. For 
example, in a hierarchical merging picture, one could instead imagine that this evolution is 
common to some but not all of the subsystems which form the Galaxy. Thus, D depletion 
(and carbon and white dwarf production) would be inhomogeneous within the protogalaxy, 
and across the observed quasar absorbers. 

Figure 1(c) shows the predicted supernova rates, both for Type II as well as Type la 
events. The model use here is that shown in the solid curve of Figure 1(a) and in Figure 
2. The Type II rate follows the assumed star formation rate, which has simply been taken 
to be ■?/' oc Mgas and which turns on sequentially after the bursting Pop HI star formation 
ends. The Type la rate is calculated using the prescription of Matteucci & Greggio (1986), 
and assumes that the binarity properties of the Pop HI stars are the same as those observed 
locally today. The Type la rate peaks sharply, echoing the rapid burst and small mass range 
that characterized the Pop HI star formation. Note that while the peaked nature of the 
la rate is a generic feature of the model we propose, the peak redshift is model-dependent, 
reflecting the characterist mass (and thus the lifetime) of the Pop HI stars. Figure 1(c) 
shows the results for a value rric = ^Mq, which gives a Type la peak at z ~ 3. Had we 
chosen rric = 6Mq {4Mq), we would find a peak at z ~ 5 (2). In any case, we predict 
that the Type la supernova rate is quite large at moderate redshift; this can be tested in 
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ongoing high-z supernova searches. Note that a suppression of Type la supernovae at low 
metalhcities (Kobayashi et al. (1998)), if real, could alter this prediction. 

A potential constraint on high-redshift Type la supernova activity also arises from the 
cosmic 7-ray background at MeV range. The dominant source of photons at these energies 
is mostly nuclear line emission due to the ^^Ni -^ ^^Co -^ ^^Fe decay chain (e.g., Ruiz- 
Lap uente, Casse, & Vangioni-Flam (2001)). The enhanced SN la production of our model 
increases the MeV background relative to the predictions expected from a constant-IMF 
cosmic star formation rate. However, most of the contribution is at high redshift and is 
thereby diluted. In fact, the background is not a good discriminant of these epochs, and the 
Type la rate proposed here is non inconsistent with the cosmic MeV background. 



4. Conclusions 

The present work is motivated by the reports of a population of halo white dwarfs 
in our Galaxy, and of dispersion in the high redshift D/H values. In light of these recent 
observations, we have explored the consequences of a high redshift population of intermediate 
mass stars. This population leads to D destruction at low metallicity, as well as white dwarf 
production, and would add to the far-IR background without increasing the local near-IR 
background. For this scenario to remain viable, it is necessary that the C and N yields of 
zero- metallicity stars are small, as suggested by some recent stellar evolution models. We 
have constructed relatively simple models which fit the data reasonably well for primordial 
abundances of D/H = 4 and 5 xlO^^. Other values are possible by merely adjusting the 
astration constant. Mi, in the POP HI star formation rate. Thus, we conclude that it is 
difficult to use existing data to accurately predict either the primordial abundance of D/H 
or n^/i^. 

An interesting prediction of a halo in which the white dwarf number exceeds that in the 
disk by an order of magnitude is that, binarity fractions being comparable, the number of 
Type la would be greatly boosted in the very early universe. Some may even be seen in the 
halos of nearby galaxies, (Smecker & Wyse 1991), though the spatial distribution of these 
events within galaxies will depend on the interplay between the formation and evolution of 
the Pop HI stars and protogalaxies. 

A potentially disturbing complication is raised by the work of Hansen (2001), who 
estimates the white dwarf ages from their luminosity and color distribution. He finds that 
they are relatively young and if true, this is inconsistent with our proposal. On the other 
hand, if the white dwarfs are real and in the halo, a young population is even more difficult 
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to explain than an old one. Given this, and the difficulty of the age calculations (which also 
couple to the difficulty in determining yields accurately), we require that the Pop III white 
dwarf relics are old. However, further study of this issue is clearly warranted. 

The recently reported observations of an intrinsic dispersion in high-redshift D, taken 
together with the observation of a new Galactic white dwarf population, demand a non- 
standard history of star formation and chemical evolution. If these observations persist, 
then an unconventional scenario of the kind we propose will be required. We note that 
standard scenarios of chemical evolution are unable to lead to significant dispersion at low 
metallicity. If these standard models prevail, then one must concur with Pettini & Bowen 
(2001) that the D dispersion at high redshift is an artifact, and that significant systematic 
errors are present. At present, the use of D/H as a precision baryometer must be called 
into question until these systematic errors are better understood. We expect that future 
analyses and observations will either show a convergence towards a unique observed D/H 
value, representative of the primordial one, or reveal the dependence of D/H on metallicity 
in high z absorbing clouds. 
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